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Abstract 

We discuss how the inckision of singular gauge fields in the partition function for 
QCD can lead to a domain-like picture for the QCD vacuum by virtue of specific 
conditions on quantum fluctuations at the singularities. With a simplifled model of 
hyperspherical domain regions with interiors of constant held strength we calculate 
the basic parameters of the QCD vacuum, the gluon condensate, topological sus- 
ceptibility, string constant and quark condensate, and briefly discuss confinement 
of dynamical quarks and gluons. 
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1 Introduction 



Singular gauge fields are most probably unavoidable in nonabelian gauge theories In 
the last two decades great effort has gone into clarifying whether or not (and, if yes, then 
how) such fields provide for confinement and chiral symmetry breaking. 

We address this question in a model based on an assumption that singularities in 
vector potentials are concentrated on hypersurfaces dVj in Euclidean space, in the vicinity 
of which gauge fields can be divided into a sum of a singular pure gauge S'^ and a regular 
fluctuation part Q^, and where a certain colour vector n'^ can be associated with dVj. 
For such configurations not to have infinite action |0 the fluctuation fields charged with 
respect to rij must obey specific boundary conditions on dVj. The interiors of these 
regions thus constitute "domains" Vj. Gauge field modes neutral with respect to n| are 
not restricted and provide for interactions between domains. In a given domain Vj the 
effect of fluctuations in neighbouring regions can be seen as an external gauge field Bj^ 
neutral with respect to n|. This enables an approximation in which different domains 
are assumed to be decoupled from each other but, to compensate, a certain mean field is 
introduced in domain interiors. 

Decomposing a general gauge field Aj^ = Sjj^ + Qj^ and demanding finiteness of the 
classical action we come to the conditions 

n^.QO^) = 0, 4j = -i fe'^'^^'ij, i) = i)i fe-''^^^\ x E dVj, (1) 

with the adjoint matrix nj = T"-n'j appearing in the condition for gluons, and a bag-like 
boundary condition arising for quarks, with a unit vector rij^{x) = x^/\x\ normal to dVj. 

To make the model analytically tractable we consider spherical domains with fixed 
radius R and approximate the mean field in Vj by a covariantly constant (anti-) self-dual 
configuration with the field strength 

= n(^)i?Ji), = ±Bl^J, Bj^jBl^ = B\„ B = const, 
nO) = cosO + sinO (0 G {^(2fc + 1), = 0, . . . , 5}), 

where the parameter B is the same for all domains and the constant matrix n'^t'^ belongs 
to the Cartan subalgebra. It should be stressed that there is no source for this field on 
the boundary and therefore it should be treated as strictly homogeneous in all further 
calculations. The homogeneity itself appears here just as an approximation. 
Thus the partition function we will deal with can be written 
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where the thermodynamic hmit assumes V, N oo but with the density = N/V taken 
fixed and finite. The measure of integration over parameters characterising domains is 
defined as 

j dai. . . = J dai J d(pi J dOi sin Oi 

3 4 5 

X r duoi ^(^^ - r £ " + • • • • 

•'^ fc=0,l ° (=0,1,2 

Here (^j and Qi are spherical angles of the chromomagnetic field, uJi is the angle between 
the chromomagnetic and chromoelectric fields, is the angle in the colour matrix n,, a, 
is the chiral angle and Zi is the centre of the domain V^. 

This partition function describes a statistical system of density N jV composed of 
noninteracting extended domain-like structures, each of which is characterised by a set 
of internal parameters and whose internal dynamics are represented by the fiuctuation 
fields. 



2 Mean Field Correlators 

In this model the connected n— point correlator 



(5°^ {X^\ 5"" ix )\ = R"t"^---"" " {X^ 

fai...a^ _ „(i)ano{i) n{j) 

''^llUl,...,^lnl^„ J "Wj"- ■ • ■ -^fllUl ■ ■ ■ -^fJ-nUn' 



of field strength tensors. 



N 



can be calculated explicitly using the measure, Eq. (Q). Translation- invariant functions 
Enixi, ...,xn) = - ! d^zOil - (xi - zf / R") ... 0(1 - (x„ - zf / R") 

V J 

emerge and can be seen as the volume of the region of overlap of n hyperspheres of 
radius R and centres (xi, . . . normalised to the volume of a single hypersphere v = 
7r^i?^/2. The functions S„ are continuous and vanish if |xj — Xj\ > 2R. Correlations 
in the background field have finite range 2R. The Fourier transform of S„ is then an 
entire analytical function and thus correlations do not have a particle interpretation. The 
statistical ensemble of background fields is not Gaussian since all connected correlators 
are independent of each other and cannot be reduced to the two-point correlations. 
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The simplest application of the above correlators gives a gluon condensate density 
which to this approximation is 

Another vacuum parameter which plays a significant role in the resolution of the Ua{^) 
problem is the topological susceptibility . To define this we consider first the topological 
charge density for the colour group SU{3) in the lowest approximation 

where uj G {0, vr} depends on the duality of the j-th domain. For a given field configura- 
tion the topological charge is additive 

Q = [ d'^xQ(x) = q(N+ - A^_), q = B'^R^/IQ, -Nq < Q < Nq 
Jv 

where g is a 'unit' topological charge, namely the absolute value of the topological charge 
of a single domain, and is the number of domains with (anti-) self-dual field, 

= A^++A^_. The probability of finding the topological charge Q in a given configuration 
is defined by the distribution 

V (O) = = ^ 

^fN 2^{N/2-Q/2q)\{N/2 + Q/2q)V 

where MNiQ) is the number of configurations with a given charge and Mn is the total 
number of configurations. The distribution is symmetric about (5 = 0, where it has a 
maximum for even. For A^ odd the maximum is at Q = ±5'- Averaged topological 
charge is zero. 

The topological susceptibility x is determined by the two-point correlator of topological 
charge density, which in the lowest approximation gives 

X = j d'x{Q{x)Qm = ^ / d'xE.ix) = (3) 

3 Area Law for the Wilson Loop 

To zeroth order in fiuctuations the Wilson loop is given by the integral 

^ f d'^z- f 1 ( f 

= v^,fcoo n i V y da,-TT exp |z J^^ da,^{x)B,,{x 



Note that path ordering in our case is not necessary since the matrices are assumed 
to be in the Cartan subalgebra. Computationally it is convenient to consider a circular 
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contour in the [x^, X4) plane of radius L with centre at the origin. To illustrate the steps 
in the calculation we consider here the case of colour group SU (2), though the final result 
for SU{3) will be quoted below. For colour SU{2) we have = t^r^^ = ±1. The 
thermodynamic limit assumes that the density N/V is fixed. Calculation of the colour 
trace gives 



W(L) = lim 



Sl 



da^,{x)9{l-{x-Zky/R''). 



N 



We have exploited the property that the integral above does not depend on the index j. 
As the contour of the Wilson loop is in the (^3, X4)-plane, the only nonzero components 



of Jfj^i, are 



J34 = -J43 = / dx3dx49{l - {x - zf/R^), Bf^^Jf^^ = 2SJ43COS6', 

where 6 is the angle between chromoelectric field and the third coordinate axis. After 
integrating over the spatial orientations of the vacuum field the Wilson loop takes the 
form 

N 



W{L) 



lim 

N,V->-oo 



— I dz 



sm2BJ43{z 



V Jv IBJ^ziz) 



Integrating over z and then taking the thermodynamic limit {N 00, v — V/N 
7r^it!^/2), gives finally for a large Wilson loop R 



W{L) 



lim 



N 



U{L) 



N 



which displays an area law. For SU{2>) the only difference is in the function f[z) which 
turns out to be: 



3772; 



3- 



V3 t'^-^^l^dx 



l-KZ 



Jo X 



smx — 



2^/3 r^/^dx \ 
/ — sm a; 
Jo X J 



TTZ 



The function / is positive for 2; > and has a maximum for z — 1.55. We choose this 
maximum to estimate the model parameters by fitting the string constant to the lattice 
result, 



^/B = 947MeV, R'^ = 760MeV, 



(4) 
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and get for the gluonic parameters of the vacuum 

= 420MeV, X = (197MeV)^ — (F^) = 0.081(GeV)^ q = 0.15. (5) 

If there were clear separation of the two scales characterising the system (namely if either 
\/BR ^ 1 or \/BR -C 1) then an approximate treatment of Green's functions for the 
fluctuation fields would be possible. We observe that there is no separation: \fBR ^ 1. 
Thus an approximation based on large or small domains is not justifiable. 

Observe also that if B goes to zero then the string constant vanishes. This underscores 
the role of the gluon condensate in the confinement of static charges. On the other hand, 
if the number of domains is fixed and the thermodynamic limit is defined as V, — >■ 
cx), = const. < oo, namely if the domains are macroscopically large, then W{L) = 1, 
which indicates the absence of a linear potential between infinitely heavy charges in a 
purely homogeneous field. 

Since we have exactly integrated over background fields the role of finite range of 
correlation functions is hidden in the above calculation. In order to see this role explicitly 
one would need to decompose the integrand into an infinite series and integrate term by 
term. At this step all correlation functions of the background field up to infinite order 
would be manifest. The arguments about the importance of a fast decay of correlators 
for confinement of static charges 0] would be seen to apply here via this representation. 

4 Confinement of Fluctuation Fields 

Confinement of the fiuctuational fields can be studied via the analytical properties of 
their propagators. For the types of fields and boundary conditions we have described 
here, the propagators of fluctuations can be analytically calculated by reduction of the 
problems for both quarks and gluons to that of a charged scalar field. The scalar case 
is essentially equivalent to the problem of a four- dimensional harmonic oscillator with 
the orbital momentum coupled to the external field, and the general solution for these 
Green's functions can be found exactly by standard decomposition over hyperspherical 
harmonics. 

Qualitatively it is clear that due to the boundary conditions we are dealing with, the 
x— space propagators of charged fields are defined in regions of finite support where they 
have integrable singularities so that their Fourier transforms are entire functions in the 
complex momentum plane. This can be treated as confinement of dynamical fields]^]. 
A known consequence of entire propagators is a Regge spectrum of relativistic bound 
states|P], which is physically appealing. A more expansive analysis of confinement in this 
model will be given elsewhere 0. 
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5 Quark Condensate 



Due to averaging over self- and anti-self-dual configurations and all possible values of the 
angle a in the partition function chiral symmetry is not broken explicitly. However, as we 
show below, a nonzero quark condensate arises in the massless limit due to an interplay 
between the random distribution of domains with self- and anti-self-dual field and the 
boundary conditions with random value of the chirality violating angle a. 
The quark propagator is defined by the equations 

(i ^ - ^fi7^-B^^x^ - m)S{x, y) = -5{x, y), (6) 

tjjix)e'''''^'Six,y) = -Six,y), {x - zf = R\ 
S{x,y)tfj{y)e-'"'''^ = Six,y), {y - zf = R\ 

where 77^(x) — {x — z)^/\x — z\. Substitution 

S{x, y)^{ip + m)[P±no + P^O+Hi + P^O^H-il 

7V± = ^(1 ± n/\n\), E± = ^(1 ± ^B/B), B = \n\B, 

shows that scalar functions T-C(;{x,y) should solve the equations: 

{-D^ + m^ + 2CB)nc = S{x,y), 

If we were to look for solutions vanishing at infinity then the Green's function Ti-i (C = 
— 1) would be divergent in the massless limit due to zero modes of the Dirac operator. 
The bag-like boundary conditions remove zero eigenvalues from the spectrum, and the 
massless limit is regular. 

In order to avoid cumbersome calculations and render the role of the former zero modes 
transparent, we turn to the particular choice y — z — and calculate the value of the 
quark condensate at the centre of the domain. In this case are functions of x^ only, 
and the general solutions for scalar Green's functions take the form (//^ = /2B -\- C) 

He = A(xVc) + Cc$(xVc)> *(^» = e-^"'/^M(l + fx, 2, SxV2). 

Here A(a;^|/i) is scalar propagator which vanishes at infinity with mass 2Bfj,, and 
is a solution to the homogeneous equation regular at = 0, expressed in terms of 
the confluent hypergeometric function. The constants can be fit to implement the 
boundary condition. The terms rnTio ^-nd mTli vanish in the massless limit and do 
not contribute to the condensate. The nontrivial contribution comes from mTl-i. The 
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bag-like conditions imply that on the boundary 7i_i satisfies a mixed condition, where 
/' = df /d\x\ and the sign (— )+ corresponds to an (anti-) self-dual domain, 



which leads to the relations 



l^^o^^-i = ^;;im^(^^V2)e-^^- /^ Tr5(0,0) = ^-^E^(^^V2), 

\n\ 

^2 .00 ^^g2t-2;(cotht-l)/2 

F(^) = e^-z-l + — / 2 (cotht-1). 

4 JO sinh t 

Note that the term in the propagator with nonzero trace is proportional to the zero mode 
of the Dirac operator: pP^O^ exp(— 5x^/4) = 0. 

Averaging this result over the angle a and (anti-) self-dual configurations and taking 
into account the a-dependence of quark determinant detS""^ oc exp{±iga}, with q 
being unit topological charge, we get a finite result 

For B and R fixed by the string constant as in Eq. (^) this is equal to 

(^^) = -(228MeV)^ 

which indicates spontaneous breaking of chiral symmetry. Final conclusions require com- 
plete calculation including averaging over positions of domains. 
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